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Several aspects of the time-dependent Schrodinger equation are discussed in 
the context of Quantum Information Theory. 



1 Time Dependent Schrodinger Equation and Quantum 
Computation 

In Quantum Information Theory (QIT), the time dependent Schrodinger 
equation pP acquires an informational content in addition to its physical con- 
tent that dictates the time evolution of a quantum system. The functioning of 
a Quantum Computer (QC) ||2| parallels the unitary evolution of a quantum 
system. Each step of the quantum evolution corresponds to a typical operation 
on a QC: 

i/ A QC starts at some initial state 

mu))--=\%)- (1) 

The preparation of the initial state corresponds to the input data in a com- 
puter program. The state \'l^{ti)) belongs to a Hilbert space of many-particle 
states 7i®" that is a tensor product of n single-particle states. When each 
single Hilbert space Ti. is spanned by a two-dimensional basis, the information 
is quantumly encoded in a qubit. A quantum register is a tensor product of 
qubit states. Quantum states can be used to encode information. A manage- 
able configuration for the initial state is a factorized tensor product of qubits 
(unentangled state) 

|f (ti)> = \xi) ® \X2) ® . . . ® \Xn) e H®", X, = 0, 1 Vl. (2) 

ii/ A QC evolves in time during an interval t allowing us to process the 
information previously encoded 



(3) 
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where the unitary evolution operator U{t,ti) satisfies the time-dependent 
Schrodinger equation, either in differential form 

ih^U{t,U)^H{t)U{t,U), (4) 
at 

or as an equivalent integral equation 
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U{t,t{) = I-- I H{t')U{t',U)dt'. (5) 



The integration of the evolution equation yields 

~jr r* H(t')dt' 

U{t,ti) ^ "J*^ 



(6) 



where T is the temporal ordering operator introduced by Dyson. When the 
Hamiltonian is time independent as in conservative systems, the integrated 
form of the evolution operator can be further simplified as 

[/(t,ti) =e-^(*-*')^. (7) 

The evolution operator of Quantum Mechanics corresponds to a quantum logic 
gate in a QC. Thus, a quantum logic gate is a unitary operator acting on the 
states of a certain set of qubits (quantum register). It is represented by a 2" x 
2" matrix in the unitary group U(2"). Hence, quantum computation is a class 
of reversible computation. The first experimental proposal for constructing 
quantum logic gates with ion traps was by Cirac and ZoUer 13, 0]- 
iii/ Quantum programming corresponds to harnessing the reversibility in- 
herent to quantum computation with the purpose of producing a useful 
task. A task corresponds to the evaluation of a certain Boolean function 
/ : {0,1}™ {Ojl}- Reversibility requires to split the quantum register 
storing an initial state into two parts: the source register and the target 
register, namely, 

l^'o) := \^s) \^t)- (8) 

To implement a Boolean function / in a QC we need the action of a unitary 
gate Uf on the registers as follows 

Uf\xiX2 . . .Xm)s\Xm+l)t = \X1X2 ■ ■ .Xm)s\Xm+l © fixi,X2, ■ ■ .,Xm))t- (9) 

This action guarantees that whichever the / is, the function evaluation is 
implemented reversibly j^. 

iv/ The final state is the outcome of a computation in a QC. It is a 
superposition of all states in a computational basis . . . , Xn)} 



^Uitf,ti)\^{ti)) = C,,,...,,Jxi,...,x„), (10) 
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and will be an entangled state in contrast to the factorized state that we 
started with Entanglement is the result of combining the superposition 
principle of Quantum Mechanics with the tensor product of single-particle 
systems. This capability of storing many classical registers (xi, . . . ,Xn) at the 
same time (|10|l is called quantum parallelism. It is the source of the potentially 
superior efficiency of a QC with respect to a classical computer, 
v/ The measurement in Quantum Mechanics is the read out of the QC to 
retrieve the desired output. Associated to the computational basis there is a 
set of orthonormal operators 7J(a;i,...,a;„) = • . • , a^n) • ■ • , a;„| such that 
the result of a measurement on the final state is 

l^^^^^^l Tr[i7(.,....,,,|^,>('^d] = ^''^ 

with probability 

7'(xi,...,a:„) = Tr[7T(,^,...,,„)|tf'f)(tf'f|] = |C,,,...,,„p. (12) 

Thus, the probabilistic nature of Quantum Mechanics implies that a quantum 
computation is probabilistic: we will find the desired output within a certain 
probability lO. The aim of a good quantum programmer is to make this 
probability as close as possible to 1. This means a pattern of constructive 
interference of amplitudes towards the desired output amplitude. Examples 
of good quantum programming are the Shor algorithm where the task is 
to factor out a large integer, and the Grover algorithm [7], where the task is 
to search for an item in a disordered database. 

This is the time evolution of a quantum system and it encompasses all 
the principles of Quantum Mechanics, which in turn can aslo be given an 
interpretation in QIT |H]. 

The basic problems that a QC can address are the same than a classical 
computer: 

• Computability: what problems the machines can do and cannot do. 

• Complexity: scaling of space/time with the size of the problems that can 
be solved. 

• Universality: whether one machine can simulate all others, so that we do 
not need to construct an special-purpose machine for each problem to be 
solved. 



2 Reversibility on the Lattice 

Reversibility means that we can compute exactly backwards by saving only 
the final conditions and the number of steps. It is energetically economical 
compared to current standard computers that work irreversibly 1^. 

A QC works reversibly as a consequence of the unitary of the evolution 
operator (0. It is also a machine that operates with finite means (Turing 
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machine): to simulate the evolution of a quantum system we must discretize 
both time and space. This is done by setting the Schrodinger equation on a 
space-time lattice. This discretization was implicit in the qubit decomposition 
of the initial state in ||2Jl. From a more fundamental point of view, many 
theoretical physicists, including Feynman have wondered about what is 
the underlying structure of the universe beyond the subatomic level and have 
proposed that it might have a discrete structure. If this were the case, a natural 
question is how to formulate reversible evolutions in a discretized form. 

The finite difference method is a natural scheme to discretize the Schrodin- 
ger equation: substitute derivatives of functions by finite differences, which 
are approximations to a certain degree in a Taylor expansion. There is some 
freedom in this transition from a differencial equation to its finite difference 
form. It turns out that some of them are not exactly reversible while some 
of them are. The meaning here of the word exactly is explained below in the 
context of computer simulations with finite precission. 

The key point is whether we use a naive difference formula for the deriva- 
tive of a given function F{z) as 

±F{z)^-[Fiz,+,)-Fiz,)] + 0{e). (13) 
az e 

This is a non-centered formula. Another possibility is to use centered formulas 
for derivatives of any degree 

^^Fiz) = l[F(z,+i) - Fiz,^,)] + 0{e% (14) 

^F(z) = ^[F(z,+i) - 2F(z,) + F{z,^,)] + ©(e^). (15) 

Let us consider the evolution of a single particle in one dimension given 
by the Schrodinger equation in coordinate representation 

^Ofix, t) = -Q^l'i^, t) + V{x, t)^{x, t), (16) 

where V{x,t) is a generic potential and we assume units H = 2m — 1. To 
discretize this equation, let a be the lattice spacing so that the discrete coor- 
dinates are Xm — ma,m € Z and r the time step such that time coordinates 
are tji = nT,n S Z. The discretized wave function is Wm,n ■= ^{xnntn) and 
similarly for the potential Vm,n '■— y(xm,tn)- 

Asymmetric Difference Scheme 

This is a (centered space, non-centered time) difference scheme that results 
from the substitution of H13() for the time derivative and (|15|l for the spatial 
derivative 
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m,n+l 



l,n ) '^^m^n^m^n 



(17) 



This is a first order difference equation in the time index n like its continuum 
version. The initial value problem is determined by giving the values ^m,o 
along the one-dimensional spatial grid. However, this equation H17|) is not 
an exactly reversible discrete update rule because of roundoff and truncation 
errors. The matrix form of H17|l is, for periodic boundary conditions in a lattice 
with M sites, 



1 



n+l 



iVi,„r - 
ie 



2ie 



ie 

'^V2,nT — 2ie ie 



le 




V 











1 - iVM,,iT - 2ie / 



(18) 

where e := ^ and ^„ is a column vector of dimension M storing the values 
of the wave function at time step n. This matrix is unitary only in the limit 
of infinitesimal e and r: the norm of each column vector is 1 + O(t^) + O(e^), 
and the scalar product of two vectors is O(t^) + O(e^). Were it be exactly 
unitary, its inversion would not need multiplication/division operations since 
it would be its Hermitian matrix: only simple arithmetic operations would do 
the job. 



Symmetric Difference Scheme 

This is a (centered space, centered time) difference scheme that results from 
the substitution of (fnji for the time derivative and lfT5|l for the spatial deriva- 
tive 



^m.n+l — ^m,n-l +i 



2^ (If ,„+!,„ - 2Wm.n + !f™-l,,i) - 2TVm,n1'm,n ■ (19) 



This scheme was discovered by Fredkin and Barton in 1975 |16| in their quest 
for exact reversibility of discretized Schrodinger equation following discussions 
with Feynman |16| . |17| . jl8) . This discretized update rule have very interesting 
properties: 

i/ It is a second order difference equation: the initial value problem is defined 
by providing the values of the wave function at two earlier time steps. This 
is in sharp contrast with the original continuous differential equation. 

ii/ It is exactly reversible after n steps of computation. 

iii/ It is exactly reversible even if the RHS cannot be computed exactly. 

iv/ It is exactly reversible even if ^rn,n has finite precision values. 

v/ It is a complex equation. However, the real part of '4'„+i depends only 
on the real part of '^n-i and the imaginary part of Likewise, the 
imaginary part of ^n+i depends only on the imaginary part of ^n-i and 
the real part of 
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The exact reversibility expressed in properties ii/,iii/ and iv/ is related to 
property i/ and follows from the backwards recurrence relation: 



m,n4-l ' 



2^ i^m+l,n ~ "^-^m^n + ^v^-l,n) ~ 2TVra,n1'ra,n ■ (20) 



Property v/ is also useful to see the exact reversibility more explicitly: we can 
select the real parts of the wave function at all even time steps n — 21 and the 
imaginary parts at all odd time steps n = 21 + 1. Thus, v/ guarantees that 
the evolution under H19|) respects this structure and the rest of components 
of *2i, "^21+1 can be ignored self-consistently. 

As a warmup, let us consider a typical update rule for a set of variables 
{X21, l2;+i}i^0' ordered according to the time evolution of the real and imag- 
inary parts described above: Xq, Yi, X2, Y3, X^, Y<^, Xq, .... The initial data is 
Xq, Yi. Let us assume that the update rule has the following form: 

X21+2 '■= X21 + F{Y2i+i), , . 

Y21+3 ■.= Y2i+i+F{X2i), ^ ' 

where F is an arbitrary function. This is an example of an exactly reversible 
recurrence relation: the backwards update rule follows inmediately, 

(Y2i+i= Y21+3 - F{X2i), 

\ X2I =X2l+2-F{Y2l+l). ^ ' 

Notice that in doing this reverse operation there is no need for multiplica- 
tion/division operations, instead the simple arithmetic of addition/substraction 
can do the job. This fact is important when implementing the exact reversibil- 
ity on a computer. 

Let us now turn to the evolution difference equations for the real and 
imaginary parts of the wave function: 

From H19|l . they are given by 



Rm,n+1 — -^m,n— 1 ['^^ {-^ni+l.n ^-^ra.n ^"-^m — l.n) ^tV^/^^ 



(24) 



For simplicity we have assumed that the potential is time-idependent Vm.n ■= 
Vm,yn. Thus, setting n — 21 + 1 and n — 21 + 2 in the first and second 
equations of (|24ll respectively, we can discard the real parts at all odd time 
steps and the imaginary parts at all even time steps, yielding 

Rm,2l+2 = Rm,2l — [2e (/m+l,2i+l " 2/„i_2i+l + Im-l,2l + l) — "^rVmlmM+l] , 
Im,2l+3 — Im,2l+1 + [2e {Rm+l,2l+2 — 2i?m,2i+2 + -Rm-l,2i+2) — '^TVmRm ,21+2] 

(25) 

They have precisely the same form as (|21ll 
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Rm,2l+2 — Rm,2l — F„i{Im,2l + l) j ^2g) 
21+3 — Im.2l + 1 + Frn{Rm,2l), 

with F.fn given by 

Therefore, we have an exactly reversible update rule for the real and imaginary 
parts of the wave function at successive time steps. This does not produce 
the values of ^„ at arbitrary times, but assuming that the change in y&n is 
smooth when r is small enough, then we can define the complex value of the 
wave function at any time as 

'l'm,l ■= Rrn.2l + ilm,2l+l- (28) 

A more compact form of the update recursion relations are obtained if we 
arrange the spatial components of the wave function in a vector array: 

'R-21 ■= {Rm,2l), ^21+1 ■= {Im,2l+l)- (29) 

Moreover, we may want to implement the evolution equations directly in 
a computer with only integer arithmetics and resorting to only summa- 
tion/substraction operations. The update rules H26|l are precisely well suited 
for this type of implementations using for instance the C/C++ languaje oper- 
ations += and -=. With these provisos, the computer-like form of the evolution 
difference equations is 

■R -= LF(i)J, 

I + = LF(R)J . ^^^> 

The total probability P{t) := ||tf/(a;, i)||2 for finding the particle is a con- 
served quantity in the evolution of the quantum system This quantity 
admits a discretized implementation like 

Pi := + l2;+i. (31) 

However, it is not exactly reversible for the discretized time evolution H26|l . 
Feynman found |16| the following invariant 

-P/''^ R-2i ~ l2i + l ■ l2i-l, (32) 

that is exactly reversible under the evolution (|26|) . 

The stability of the exactly reversible Schrodinger equation is analysed in 
appendix IXI 

So far we have considered the reversible vs. non-reversible implementation 
of the time-dependent Schrodinger equation in a classical computer. We can 
also implement it in a QC. Feynman proposed that one of the most 
important tasks that a QC could carry out is simulating quantum physics in 
an efficient way, contrary to the inefficiency showed by classical computers. 
The Abrahams-Lloyd scheme can serve for implementing differential equations 
in a QC EOl- 
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3 The Schrodinger Equation and the Grover Algorithm 

The initial steps pursued by Schrodinger to formulate his evolution equation 
for the wave function of a quantum mechanical system bears little similar- 
ity with the final outcome he arrived at. In fact, he first came up with a 
scalar relativistic wave equation, later known as the Klein-Gordon equation. 
The proccess of creation of a new theory usually is not a straight line but 
rather a long a winding road, whose final results are not as polished and neat 
as they look at the end. Likewise, the final formulation of Grover algorithm 
does not show the initial steps that his author followed. In fact, Grover wan- 
dered around the discretized time-dependent Schrodinger equation |5J until 
he arrived to his beutiful algorithm in its final form [Jj, [21], |22j . 

The Grover algorithm is a quantum searching algorithm that outperforms 
any classical searching algorithm. The problem of searching an item in a disor- 
dered list or database with M items is known to be one of the most fundamen- 
tal tasks that any computer is oftenly doing. Classically, no algorithm can do 
better than a brute force algorithm 0. The number of queries is order 0{M). 
Remarkly enough, the quantum search algorithm can do better with the help 
of the quantum parallelism and constructive interference of amplitudes in the 
evolution of a QC described in Sect. In the quantum case, the number of 
queries is reduced to order 0{y/M), which is a major improvement. 

The basic steps that led Grover to the quantum search algorithm starting 
from the discretized time-dependent Schrodinger equation are the following: 

1/ Discretize the time-dependent Schrodinger equation on a space-time lat- 
tice with lattice spacing a and time steps of r duration as described in 
Sect. 12 using an asymmetric difference scheme (|17|l . (|18f) . This evolution 
is approximately unitary up to order O(e^) + O(t^). 

It is interesting to notice that Grover does not resort to an exactly reversible 
implementation of the evolution equation as in a symmetric difference scheme 
1)19(1 . He employs a different pathway to achieve an exact unitary transfor- 
mation since he needs to stick to a first order difference equation in the time 
variable. 

2/ Introduce Diffusion D and Rotation R operators. This is a decomposition 
of the matrix evolution equation H18|) as follows 



n+l 



DR^ 



(33) 



/ 1 - 2ie ie 0.. 
ie 1 - 2ie ie . . 







D : 



(34) 







0.. 



1 - 2ie / 
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R : = 





-iV2T 



\ 





\ ... e-'^^^^ / 

We have assumed that the potential is time independent. This decompo- 
sition is correct within the approximation of order O(e^) + O(r^). Notice 
also that while R is exactly unitary, D is only within the approximation 
0{e^). 

The motivation for doing this decomposition is the key idea behind the con- 
struction of the searching algorithm: the matrix D may be interpreted as a 
diffusion transformation (Markov process) with imaginary transition proba- 
bilities ie, and the matrix i? is a phase rotation at every site of the lattice 
m = 1, 2, . . . , M induced by the application of the external potential Vm- 

In the searching process, at time step n — the system is initially in an 
equally probable superposition of states from the database to be searched for 



(35) 



1 



1,2, 



,Af. 



(36) 



The goal is to find a marked state or item that is represented by a marked 
component of the wave function ^ma,N, after a number N of iterations that 
represent the queries of the disordered database. This can be achieved taming 
these transformations D, R in order to drive the system towards a final wave 
function that has the marked component highly amplified with respect to the 
rest of components. 

The diffusion process represented by D affects only to nearest-neighbour 
components of the wave function. This can also be viewed as the hopping of 
a single particle along the lattice in the presence of a potential. A classical 
analogy is useful: we may think that this particle will roll down towards the 
sites of lower potential like in a potential well. This is the amplification effect 
needed in the quantum search. 

3/ Global diffusion and localized rotation: it is advantageous to extend the 
local diffusion of eq. H14() globally, as follows 



( 



\ 



{M-\)\e ie ie 

ie 1-{M - l)ie ie 



le 
ie 



(37) 



le 



le 



1 - (Af - l)ie / 



Now, the particle may hop to any site of the lattice. This is an example 
of long-range hopping. The sum of every column in (|37|l is 1. 
At the same time, we tune the form of the potential well to be localized 
at the marked state as follows 
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/I 

e-'"^ 
1 



VO 








(38) 



where v is the vahie of the potential at that site. When it is negative it 
means attraction. The marked state is unknown, but we can selectively 
rotate its phase by an amount given by v. 

With a global diffusion operation H37() there is no net hopping between 
unmarked sites since they have the same amplitude and the hopping in both 
directions cancels out. The hopping occurs towards the marked site from the 
rest of all states. This is the mechanism of amplitude amplification in the 
search algorithm. This amplification is maximum for a value of the rotation 
parameter v ^ ^. After this rotation and at time step n, let the amplitudes 
of the marked and unmarked states be given by 



M 



VM' 
= 1. 



/M 



(39) 



Initially, Co 1 cq (|36|I . After the diffusion transformation, they become 



1'n 



(40) 



This means that the marked state changes its magnitude an amount of the 
order 0{y/Me), since by normalization i < c„ < 1, assuming C„ < \J^- 
We wish to make e as large as possible in order to get the highest magnitude 
change in the marked state. However, we recall that the diffusion operation 
H37|) is unitary assuming that e ^ jj. 

As the change of the marked state is order 0{^/ Me), if we take e = 0{jj), 
then this change becomes 0{M~i). Thus, after the application of the diffusion 
transformation H37|l a number of time steps N = 0{\/M), the amplification 
of the marked state will become of the order of 1. 

4/ Exactly unitary diffusion operation. Let us get rid off the infinitesimal 
character of the diffusion 1)37(1 introducing a finite Dp transformation 



f X y y y .. . y\ 
y xyy ... y 

\y y y y . . . X J 



(41) 



where x, y are complex amplitudes satisfying the unitary constraints 
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|x|^ + (M-l)|y|^:=l, {xy* +x*y) + {M -2)\y\' -.^Q. (42) 

The infinitesimal transformation (|37|) corresponds to the particular case 
X = \ — {M — l)ie, y ~ ie. Solving (|42|) we find that the highest value of 
y can be obtained on real solutions as 

2 2 

X = -1 + ■■= XM, y=j^--^VM- (43) 

Then, the best diffusion operation is Df(xm jVm) which corresponds to a 
hopping amplitude as high as 

This is the pathway followed by Grover in order to implement an exactly 
reversible evolution while retaining the first order character of the iteration 
process. 

5/ The quantum sarching algorithm: iterate a number N of time steps the 
basic operations (|41() and H38() . starting from the uniformily distributed 
state (|36ll . When the diffusion operation is exactly unitary (|41|l . 142|) . the 
maximum hopping into the marked site occurs for a rotation of v — t: 
(I38|l . instead oi v = ^ that is the case when the diffusion is infinitesimal 
(123 . Thus, 

= [DF{xM,yM)RLMf *o- (44) 

This is an example of Trotter-Suzuki transformation j23| . |24j . The count- 
ing of the number of iterations -queries of the database- is the same as 
in iv/ but with the application of H41|) . H42|l to the state 



^2 



+ = 1 



The result is, 



M 



'Z^mo,n - ^ {Cn + - l)c„ ^- C„ 



(45) 



(46) 



In the limit M — > cxd, the increase in the amplitude of the marked state 
is of order 0{^^), and the amplitude of the unmarked states remains 
approximately the same. As c„ remains order 0(1) during the evolution, 
the net amplification is order 0(M~2). Therefore, we need a number of 
iterations (|44|l of the order N = 0{y/M) to make the probability of success 
for finding the searched item close to 1. 

The origin of the higher efficiency of the quantum search algorithm over 
the classical algorithms is thus a neat combination of quantum paralellism 
and constructive interference of amplitudes. In addition, the Grover algorithm 
H44|) is an example of the evolution of a quantum system like in Sect. ^ with 
^^AT — ^ if'f and SI/o — > in order to implement a useful task in Quantum 
Information Theory. 
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4 Conclusions 

It has been two decades now since the BB84 quantum cryptography protocol 
was presented P3 ^^'^ decade ago that Shor algorithm was introduced [HI ■ 
I has been specially during this last decade that Quantum Information Theory 
has emerged an matured as a solid and complete discipline. QIT has part of 
its roots in the fundamentals of Quantum Mechanics, and thus it has become 
a part of Theoretical Physics by now. It has ramifications in other disciplines 
like Information Theory, Computer Science etc. which guarantees that findings 
in QIT may have applications to other fields beyond. Specially rewarding is 
the fact that QIT is a discipline in which there are experiments testing its 
theoretical findings. All these aspects have made QIT a very exciting field so 
far and we expect that it will continue to be so in the future. 

Dedicated to Prof. Alberto Galindo Tixaire on occasion of his 70"^ birthday. 
Acknowledgments. This work is partially supported by the DGES under con- 
tract BFM2003-05316-C02-01. 



A Numerical Stability of the Exactly Reversible 
Schrodinger Equation 

In the analysis of the finite difference equations obtained from the time- 
dependent Schrodinger equation Hl()|l . who have focused on the exact re- 
versibility on the lattice as the main issue. 

Another very important issue is whether the difference equations are well- 
posed: a solution of the difference equation must be a reasonable approxima- 
tion to the solution of the differential equation. This is the notion of conver- 
gence: in the limit of both the lattice spacing a — > and the time step r — > 0, 
the solution of the difference equation must converge to the solution of the 
differential equation. The convergence is measured with a L^-norm: 



A related issue is the consistency of the difference equation. It means 
that the difference equation itself must converge to the differential equation 
pointwise at each grid point as a,T 0. Consistency implies that a smooth 
solution of the differential equation is an approximate solution of the difference 
equation. While a direct proof of convergence is a rather difficult problem, 
the proof of consistency is simpler. For the exactly reversible Schrodinger 
difference equation H19(l . consistency is a direct consequence of the Taylor 
expansions H14|l . (|15|l and it is checked directly. 

Consistency is a necessary condition for convergence, but it is not sufficient. 
Additional properties must be required for a consistent difference equation to 
be convergent. This is where the notion of numerical stability enters JHl, 




(47) 
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The essence of stability is that the solutions of the difference equation must 
be bounded is some sense, in order to prevent the finite solution to explode in 
the continuum limit. More precisely, a difference equation is stable when the 
solution at a time step n is bounded by 

\\^n\\l<CTY.\\^n'\\l (48) 
n'=0 

with Ct some constant, and iVo some integer that depends on the order of the 
difference equation. The meaning of this bound is that as time goes by, the 
growth of the solution is limited. This concept of stability is closely related to 
the concept of well-posedness of initial value problems for partial differential 
equations: 

/+00 r + co 

\<P{x,t)\^dx <Ct \1'{0,t)\^dx, te[0,T]. (49) 
-OO J —OQ 

The fundamental theorem in the theory of finite difference equations is the 
Lax-Richtmyer equivalence theorem: a consistent finite difference equation for 
a partial differential equation, for which the initial value problem is well-posed, 
is convergent if and only if it is stable. 

This theorem is a very powerful characterization of convergence since the 
determination of consistency and stability is much easier. Let us apply the 
Lax-Richtmyer theorem to the exactly reversible Schrodinger difference equa- 
tion H19fl with vanishing external potential. We already know that the equation 
is consistent from the Taylor expansions (|14|l . H15() . Let us focus on stability. 
In this particular case, we have translational invariance and we also assume 
an infinite spatial lattice. This allows us to Fourier transform to momentum 
space the solution tfVt^„ in coordinate space, 



if„(fc):=-= > e-"""'=if™,„a, fce[— ,-]. (50) 
J2tt a a 

^ m— — oo 

Using the inversion formula 

1 f+» 

^,n,n = ^ y „ e'™'^'=tf'„(fc)dfc, (51) 

we obtain the recursion relation satisfied by the Fourier components 
#„+i(fc) = #„_i(fc) - 2i/,(fca)!f'„(fc), 



feika) :=4esin- 



2 



This is a simpler recursion relation since it is only on the time-step label n. 
The Parseval relation 
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+00 



^ |#„(fc)pdfc==^|tf',„,„pa=||*„||2, (53) 

"■1 -00 

guarantees that the stability property is preserved under the Fourier trans- 
form. There are several ways for solving 152|l . Let us introduce a generating 
function 

oc 

GfeH :=^*'„(fc)u", (54) 

n=0 

defined on the complex plane u G C. Knowing Gk{u), we can recover the 
Fourier components through the Cauchy theorem, 

4(fc) = 7^-^ Gk{u)u-"-^du, (55) 

where Fq is a contour encircling the origin counterclockwise. Using the recur- 
sion relation we can find the generating function. It is given by 

^ [l + 2iMka)u]%ik)+u<Pi{k) 
^'^("^ = l + 2iMka)u-u^ ' 

where ^^{k) ^Wi^k) are the Fourier transform of the initial data (|5U|I . To find 
the solution via (|55|l it is convenient to change variables w := i, 



Mk)^^i Gfe(-)t."+i^, (57) 



where is a contour encircling the point of infinty and enclosing all the 
poles of the generating function. Using Cauchy's residue's theorem, we can 
obtain the general solution as 

•^„(fc) = 2^ Hk{v)v-dv = EpoicsRes[Hfe(t;)z;"], 

TT [^'+2i/,(fca)]-|-o(fc)+^i(fc) loSj 

^kK^J) •— v'''+2if,{ka)v-l ' 

This solution depends on the roots of the equation 

+ 2i/,(fca)w -1 = 0, 



v±=-iMka)±^l-f^ika). 



(59) 



When the roots are distinct u+ 7^ W-, we have simple poles in ()58|l and the 
solution has the following form 

Mk) = \ ' + 2^Uka))Mk) + Aik)] <. (60) 



2^1 -/.^M 

The roots are equal v+ := vq —: v^, when feika) = ±1 and then vq = ^i. 
Thus, we have a double pole in H58|l and the solution takes the form 
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'Pnik) = \{n - l)#o(fc) T m#i(fc)j e=F'("-2)^. (61) 

The stability analysis can be done with the help of the explicit solutions. 
To simplify, let us assume that the limit a, r — > is taken with e kept fixed 
|26|. When there are two different roots, the solution H()U|I is bounded when 

\v±ika)\ < 1, (62) 

and this determines the stability. If 4|e| > 1, then there exists values of ka 
for which |f±(fca)| > 1 and becomes unstable. It is when 4|e| < 1 that we can 
guarantee condition (|62|l Vfca, since then |i;+(fca)| = |w_(fca)| = 1. 

When there is a single root, the solution (|61|l is not bounded due to the 
linear growth with n. This instability occurs at 4|e| sin^ ^ = 1. Since from 
the previous case we already know that 4|e| can be at most 1, removing this 
possibility we avoid the case of solution (|61|l . Therefore, the condition for 
numerical stability is 4|e| < 1. In this case, the solution (|60|l can be further 
simplified as 

'^"(^) = o .Ifu ^ E W^'-^'^^'^Mk) + sAik)] s«e'(i-^)"^'('='^), (63) 

^ cos t/g I f£Ci 1 L 

where we have introduced the following parametrization 

/e(fca) sine',(fca), 4|e| < 1. (64) 

Let us point out that the stability analysis carried out here can be repeated 
for the exactly reversible Schrodinger equation in terms of real and imaginary 
parts H26|l . H27|l . In this case we would need to introduce a couple of generating 
functions, one for each component. 
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such that f{x\, . . . ,Xm) ~ f{yi,...,ym)- Choose the states as orthonormal 
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8. The Principles of Quantum Mechanics pP| admit an interpretation in the hght 
of Quantum Information Theory as follows: 

• PI States: To each physical system there corresponds a complex and sepa- 
rable Hilbert space Ti of states. A pure state at instant t is represented by a 
unitary ray \<I'{t))B. € Ti.. A mixed state is represented by a density matrix 
operator p : = p, Tip — 1, p >0. 

In QIT this corresponds to encoding information in a quantum way in 
those states. The superposition principle encoded in the Hilbert space struc- 
ture implies the quantum paralellism for the information stored. When the 
Hilbert space Ti, is constructed by tensor product of single-particle states, 
there appears entanglement that allows for non-classical capabilities of in- 
formation processing and communication. 

• PII Observables: Each observable of a physical system is represented by a 
linear self-adjoint operator A acting on the Hilbert space of states Ti.: 

{iflA^p) = (A^'ipltP), Wtl^ e D{A),y^ e D(A''). (66) 

In QIT this corresponds to the particular physical realization of a compu- 
tational basis associated to the spectral decomposition of the operator A. 
For instance, it may be the third component of the spin of a spin-i particle. 

• PHI Measurements: If a physical system is in a pure state \^) G Ti, normal- 
ized, then the probability of obtaining a value a when measuring a certain 
observable A is 

r^{a)^\\E^mf, (67) 

where {Ea} is an orthonomal set of projectors El — Ea : EaEt — SatEa, 
Ea ~ I, associated to the spectral decomposition of A : A\a} — a\a), 
Ea ~ Generalized measurements are described with POVMs. 

In QIT this corresponds to the probabilistic nature for the result of a quan- 
tum calculation. It implies that good quantum programming amounts to 
maximizing that probability above a finite threshold. This is summarized 
as the need for constructive interference of final amplitudes towards the 
desired output amplitude. 

• PIV Collapse: If a given physical system is in a state p, the result state of 
an ideal measurement of an observable A, with respect to a set of values in 
a set A, is described by the following density matrix 

P^P-'^^^ Tr[pL{A)] T.^-P''- (68) 
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In QIT this corresponds to the final state of the QC after the computation 
has finished. 

• PV Schrodinger Equation: Between two consecutives measurements, the 
pure states of a physical system remain pures, and there exists in each 
unit ray |<f'(t)R) some representative state vector \^(t)) such that the time 
evolution is given by the Schrodinger Equation 

m^^Yf{t))=H{t)W{t)). (69) 

In QIT this corresponds to quantum logic gates represented as unitary evo- 
lution operators and the reversibility of quantum computation. The condi- 
tional logic needed to process information is achieved by acting on Hilbert 
spaces with at least two-particle states. This conditional logic is realized 
in the Hamiltonian H{t) of the system, as the physical interaction between 
two parts or particles in the Hilbert space of states. 

• PVI Quantization Rules: In a physical system with cartesian coordinates 
of position qi,q2, . ■ . ,qN and conjugate momenta pi,P2, ■ ■ ■ ,Pn, then the 
corresponding operators Xr, Pr that represent those observables must fulfill 
the commutation rules in the Weyl form 

UcVp = e-'^VffUc. Va„ A G R (70) 

[/„-e-*"-^ \/^-e-^^P 

If the system has an observable whose classical expression is A{qr,Pr,t), 
then the corresponding operator is obtained from this expression, written 
in a convenient form, upon substitution of the variables qr,Ps by Xr, Ps- 
In QIT this corresponds to the transition from a classical computation to a 
quantum computation. A phyical system can be used to compute classically 
or quantumly depending of its physical description. This is the meaning of 
information is physical. 
9. Irreversibility implies heat dissipation according to the second law of thermody- 
namics: information cannot be lost and every erased or lost bit of information 
is shed as a tiny pufi^ of heat. This shedding of information is happening perma- 
nently at each step of the fast and long chains of computations. Intermediate 
results in these calculations are discarded since they are apparently useless and 
the computer needs to reset the memory registers. This erasure produces heat 
that must be removed from the machine to prevent from malfunctioning. The 
nature of this heating is fundamental, intrinsic to information erasure. 
There are other sources of heating in a computer like the dragging of electrons 
in the currents inside the wires due to resistance. In actual computers this heat 
dissipation due to resistance is larger by far than the heating due to information 
erasure. There is still a lot of room for reducing resistance heating, for instance, 
introducing superconducting electron currents. When this would be achieved, 
then we face the heating due to dumping information. Landauer first quantified 
|1U| in energy this loss per bit as 

AE = kBT\n2. (71) 

This is an irreducible amount of heat that must unavoidably be shed when a 
computer clears its registers by erasing a bit of information. The less infor- 
mation systems erase, the less power they use. Bennett found a way to 
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overcome this heating introducing the concept of a reversible computer that 
is universal as a Turing machine. Lecerf had introduced |12| this idea of a 
reversible Turing machine earlier in the sixties. In the Bennett machine, infor- 
mation is never throw out, thus causing no heat. The problem then becomes 
one of a machine with an infinitely large memory. To continue working, the 
computer needs somehow to be reset and here is where reversibility comes in: 
rewind the machine till the initial state without erasing any single bit. These 
computers operate with reversible algorithms that do not require data to be 
erased. Bennett designed a memory-saving method in which a computer would 
carry out a few steps of the calculation, copy the result and rewind. This is at 
the expense of processing time. 

There has been recent interest in reversible classical computers for reducing a 
system's heat level while allowing chip performance to scale. The first reversible 
computer is a chip called Flattop |13j developed at the MIT's Reversible Com- 
puting research group. It is based on the reversible and adiabatic transfer of 
charge in digital circuits implemented in standard silicon. Flattop implements 
the Billiard Ball Cellular Automaton 14 , a univeral and reversible model of 
computation introduced by Fredkin. It is not a practical architecture for per- 
forming arbitrary computations, but it is regarded as a proof-of-concept re- 
versible machine showing that universal reversible computers may be accessible 
in the near future. 
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